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ABSTRACT 



We calculate the expected white dwarf luminosity functions and discovery 
functions in photometric passbands, if these stellar remnants provide a 
substantial fraction of the sought Galactic dark matter, as suggested on various 
observational grounds. We demonstrate the extremely rapid variation of the 
photometic signature of halo white dwarfs with time and thus the powerful 
diagnostic of white dwarf colors to determine the age of the Galactic halo. 
We also consider the various indirect constraints implied by a white dwarf 
dominated halo. These calculations will guide present and future observational 
projects at faint magnitudes. This will enable us to determine not only the 
nature of the Galactic dark matter but also the age and the initial conditions of 
the Galaxy formation. 



Subject headings: dark matter — Galaxy: halo — Galaxy: stellar content - 
white dwarfs 
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1. Introduction 

The nature of Galactic dark matter remains a major puzzle. The MACHO collaboration 
has now more than 12 events observed towards the LMC (Alcock et al., 1998). They 
confirm an optical depth r = 2. 91^9 x 10~ 7 and a mass of dark objects within 50 kpc 
M 50 = 2.0iJ;? x 10 11 M Q , at least 20% of the mass of the Galaxy within this limit. The time 
distribution of the microlensing events is narrow ranged (t ~ 20-70 days) with (t e ) « 40 
days (for the Einstein radius crossing time). Assuming a standard 1/r 2 isothermal halo, 
this means that a fraction ~ 60 ± 20% of the sought Galactic dark mass might be in the 
form of objects with an average mass (to) = 0.5ioi2 M & . 

Brown dwarfs and low-mass stars (m<lM ) are excluded as a significant halo 
population both by the microlensing time distribution and by the star count analysis of 
the Hubble Deep Field (HDF) at faint magnitudes, yielding for the dark halo a maximum 
stellar plus brown dwarf density p^ +SD < 0.001 x p h , where p h « 10~ 2 M Q pc~ 3 is the dark 
halo local dynamical density (Chabrier and Mera, 1997; Gould, Flynn and Bahcall, 1998). 
It is important to stress that both these observational constraints on the stellar content of 
the Galactic dark halo and the narrow range of time distribution of the events observed 
towards the LMC imply an initial mass function (IMF) different from a Salpeter one below 
~ 1 M . The recent determination of the Galactic local volume density (Creze et al, 
1998) leaves essentially no room for dark matter in the disk, and thus strongly favors a 
non-dissipative component for the Galactic dark matter. The suggestion that the events 
might be due to an intervening stellar population (dwarf galaxy or tidal debris) along 
the light-of-sight (Zaritsky and Lin, 1997; Zhao, 1998) is still controversial (Beaulieu and 
Sackett, 1998; Alcock et al., 1997). Under these circumstances, WD's remain the most 
favorable candidates for the observed microlensing events and the Galactic baryonic dark 
matter, in particular after the recent demonstration that some of the faint blue objects in the 
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HDF are consistent with very cool (T eS <4000 K) H-rich atmosphere WD's (Hansen, 1998). 
In this paper we calculate the expected halo white dwarf luminosity functions (WDLF) and 
discovery functions for different halo ages, in various photometric passbands, for a halo WD 
population consistent with the dark mass inferred from the MACHO observations. This 
works extends beyond our initial work (Chabrier et al., 1996) by including new atmosphere 
models and spectral colors appropriate for the halo population. 

2. Halo white dwarf luminosity function 

For WD's to provide a mass fraction Xwd of the halo local dynamical density ph, the 
WDLF must be normalized as : 

jndM bol = ^^xp h (1) 
The WDLF reads (Chabrier et al., 1996; Adams & Laughlin, 1996): 

dn rm sup(L ) dt coo i(m,L) dt coo i dm 

GfMbol Jm mf(L) oM ho i oM bo i dt 

where m is the WD progenitor mass and if)(t) and (f>(m) denote the stellar formation rate 
(SFR) and the IMF. The second equality stems from the assumption that the star formation 
burst at the early epoch of the Galaxy is much shorter than the age of the halo, so that the 
SFR can be approximated by a Dirac function ip(t) = S(t — t ). In that case v(t) is simply 
the number of stars such that tMs + tcooi = hi-, where t h is the halo age and t = t MS + t coo i 
is the WD total age, i.e. its cooling time plus the main sequence lifetime of its progenitor. 

The IMF in (2) must fulfill several constraints: the HDF observations imply 
P(m<i.o) << 0.01 x p h (see above) and the presence of type II supernovae at finite redshift 
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(Miralda-Escude and Rees, 1997) imply a finite fraction of stars above ~ 8M . We elected 
a truncated power-law function (Larson, 1986; Chabrier et al., 1996): 

Mm) = 4— = Aexp'^ 13 m~ a (3) 
dm 

Equation (3) approaches a power-law form m~ a at large masses and can easily be 
ajusted to reproduce any observed SNII rate, while the integral of the mass function, 
which is what matters in the present context, is determined essentially by the peak 
m p = (P/a) 1 ^ m and by the parameter (3. In order to examine the dependence of the results 
upon the IMF parameters, the present calculations have been conducted with m = 3.5 M , 
(3 = 3.0, a = 5.0, hereafter IMF1, which yields an average WD mass (mwD) ~ 0.8 M© and 
m = 2.4 M Q , (3 = 3.0, a = 5.0, hereafter IMF2, which yields (m WD ) w 0.7 M . The slightly 
larger average mass for halo WD's than for disk WD's is motivated by (i) the smaller 
mass-loss during evolution for metal-poor stars (Maeder, 1992) and (ii) by the fact that the 
faintest observed disk WD's have masses ~ 0.7-0.8 M (Leggett, Ruiz and Bergeron, 1998). 
Both IMF1 and IMF2 yield a mass-to-light ratio M/L >> 100, as required for a dominantly 
baryonic halo. Note in passing that this type of IMF provides a natural explanation for the 
lack of zero-like metallicity stars in the Galaxy, the so-called G-dwarf problem. 

The present calculations include the most recent atmosphere profiles and synthetic 
spectra calculations for pure hydrogen atmosphere (so-called DA) WD's (Saumon and 
Jacobson, 1998), and the most updated WD interior physics, C/O profiles (Salaris et al., 
1997), equation of state and crystallization along evolution (Segretain et al., 1994; Chabrier 
et al., 1996). To illustrate the rapid cooling of halo WD's, we found out that although for 
t h < 12 Gyr the entire WD population is brighter than My = 20, i.e. M bo i ~ 21-22, after 
14, 15 and 16 Gyr, only ~ 80%, 60% and 25%, respectively, of the total WD population 
remains brighter than this magnitude. The more massive WD's have cooled fainter. For 
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pure-He (more transparent) atmosphere WD's, the situation is even more dramatic and 
after ~ 8 Gyr, the majority of these stars have cooled fainter than M^ Q \ = 21 and will 
thus escape detection. However, using the Alcock and Illarionov (1980) accretion formula, 
dm/dt « 10 -20 ( o ™ M ) M & .yi~~ l , these stars will accrete ~ 1O _13 M of hydrogen, i.e. about 
a photosphere mass, during each disk crossing and will thus very likely cool like H-rich 
atmosphere WD's. 

Figure 1 displays the expected halo WDLF for DA's with the IMF1 (solid line), 
normalized to Xwd = 50%, in the most favored V-band, for th = 14, 15 and 16 Gyr, 
about the age of the oldest globular clusters (Pont et al., 1998)Q. The end of the observed 
disk WDLF is also displayed, as well as the observed WD's with halo-like kinematics, 
i.e. Vtan > 250km.s _1 (Liebert, Dahn and Monet, 1988). It is important to mention the 
sensitivity of the calculations upon the different input parameters, namely (i) the progenitor 
mass - WD mass relationship, (ii) the progenitor main sequence lifetime and (hi) the IMF. 
All the present calculations were done with the disk characteristic relationships (Iben and 
Tutukov, 1984). Although modifications of the two first relationships were found to affect 
only moderately the WDLF, the IMF is determinant. This is illustrated by the dashed line 
in Fig. 1, which displays the results of the same calculations with the IMF2. Although the 
peak of the WDLF is not affected significantly (shifted by ~0.5-l mag), the bright part 
of the WDLF, which stems from the low-mass tail of the IMF, has changed by orders of 
magnitude. This shows the extreme sensitivity of the bright part of the halo WDLF to the 
ill-constrained low-mass end of the IMF and thus the necessity to observed the bulk of the 
halo WDLF in order to constrain the IMF of the Galactic halo. 

Figure 2 shows WD evolution in a color-magnitude diagram (CMD) and illustrates 

^-Note that t~15 Gyr corresponds to the age of the universe for i7 ma4ter =0.24, 0^=0.62 
and H =70 km s -1 Mpc -1 , the presently most favored parameters. 
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the rapid variation of halo WD optical colors with time. The photometric observation 
of identified halo WD's will thus provide a powerful diagnostic to determine the age of 
the Galactic halo. It also provides complementary information to the observation of the 
WDLF. As shown in Fig. 1, a very small number of objects at My < 20 might reflect either 
a negligible fraction of halo WD's or a halo age older than 16 Gyr. The CMD will resolve 
this ambiguity by providing a determination of the age of the objects. 



3. Indirect constraints 

A Galatic halo composed dominantly of stellar remnants like WD's implies several 
constraints on the Galaxy genesis and evolution. The mass fraction of returned gas between 
the progenitor and the WD with respect to the remnant mass can be written (still assuming 
a Dirac-function for the SFR): 



J~o 8 5 M { m ~ m WD (m)}(f)(m ) dm 
J~o 8 5M m WD (f)(m)dm 



The IMF1 and IMF2 yield R(t) = 3.5 and 2.7, respectively. The total mass of 
returned gas is thus M gas 1.5 x 10 12 M Q if X WD = 0.5. Most of this gas will appear 
in the Lyman-a forest and the total mass of baryons in the leftover gas and in the 
stellar remnants is bound by the total baryon density, for the same redshift value : 
^rem(z) + Vt Lya {z) < Vt B « (0.019 ± 0.002)/i~ 2 (Fields et al., 1998). Present determinations 
yield &Ly a ~ 0.02 at z ~ 2 (Petitjean et al., 1993). The leftover gas is likely to be ejected in 
the intergalactic medium (IGM) by a wind due to the primordial generation of SNII or by 
the more efficient merging mechanism (Gnedin, 1998). This implies the presence of some 
hot (~ 10 5 -10 6 K) X-ray gas in the Local Group (LG). Although there are only hints for 
the presence of such gas in the LG, its presence has been established in other galaxy groups 
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and in the intracluster medium (see e.g. Mushotzky et al., 1996). The total mass of metals 
ejected during the envelope ejection phase of the WD progenitors is : 



where yz denotes the metal stellar yields. As noted by Gibson and Mould (1997), a strongly 
peaked IMF around 2 M will produce [C, N/ 0} abundances during the AGB phase about 
10 times the observed value in halo stars. This result, however, depends entirely on the 
assumed stellar yields yz- Gibson and Mould used yields appropriate for Z>1CT 2 x Z Q . 
Stellar evolution calculations for zero-like metallicities (Chieffi and Tornambe, 1984; 
Fujimoto et al., 1984) show that for a central degenerate core M>O.77M , thermal pulses 
along the asymptotic giant branch do not occur, so that the bottom of the convective 
envelope never reaches the carbon-enriched region and remains unpolluted. This leads to 
no C-enhancement of the interstellar medium during the planetary nebulae phase. For 
zero-like metallicities, this core mass corresponds to m>3M , a condition statisfied for 
most of the stars with IMF1 and marginally satisfied by the IMF2. The C-enrichment 
constraint might thus be relaxed for primordial stars, depending on the IMF. On the other 
hand, the presence of an unexpected level of heavy-element enrichment in the IGM at high 
redshift has been established observationaly (Cowie and Songaila, 1998). The identification 
of AGB stars as the origin of some of these elements would bring immediate support to the 
halo WD scenario. 

The halo WD scenario seemed to have been excluded on the basis of the observed 
rate of type la SN in galaxies (Smecker & Wyse, 1991). These calculations, however, 
assumed that the merging of two WD's whose total mass exceeds the Chandrasekhar 
mass would produce a SNIa. Recent calculations (Segretain, Chabrier and Mochkovitch, 
1997) seem to exclude, or at least strongly unfavor the formation of SNI by this scenario. 
These calculations seem to be supported on various observational grounds, suggesting that 
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the rate of SNIa from merging WD's has been significantly overestimated (see Segretain 
et al., 1997; Maxted and Marsh, 1998). Note also that usual assumptions about binary 
parameters in the Galactic disk (mass loss, orbital radius, rate) are likely to be irrelevant 
under the completely different primordial halo conditions. Chariot and Silk (1995) showed 
that a halo WD population such that X WD > 10% would correspond to a progenitor light 
at redshift z < 3.5 at odd with the observational constraints. These calculations, however, 
were done for a halo age th =13 Gyr and for stellar evolution models with solar metallicity. 
Low metallicity stars are brighter for the same mass and thus evolve more quickly. This - 
and an older halo age - weakens the Chariot-Silk constraint or at least pushes it to larger 
redshifts. The background light of the progenitors of a WD-dominated halo is constrained 
also by the total amount of energy distribution in the universe determined by the DIRBE 
observations (Guiderdoni et al, 1997). The IMF1 (resp. IMF2) yields (m) ~4M (resp. 
~ 3M ) for the progenitors, i.e. (L) ~ 80 x L© (resp. ~ 30 x L & ) over ~1 Gyr. Since the 
mass of the LG is M LG ~ 2 x 10 12 M , this yields (L LG ) < 10 47 erg/s. The radius of the LG 
protogalactic region was R ~ 1 Mpc (Peebles, 1993), which yields a surface brightness at 
redshift z, [i z < (2.6 x 10~ 4 )/(1 + z) 4 erg s -1 cm -2 sr _1 . For a halo formed at redshift z >4, 
this yields a peak distribution /i z <4.0 x 10~ 7 erg s _1 cm~ 2 sr~ x at \ z = A (l + z) > 1 fj,m, 
about a factor 50 below the DIRBE limit. 



The number of WD's of absolute magnitude M A per arcmin 2 observable in a field of 
longitude and latitude (/,&), for a limit magnitude mi reads: 



4. 



Predicted counts 



N W D = 
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v 180 ; V 
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where R = 8.5 kpc is the galactocentric position of the Sun and log d max (M\) [pc] = 
0.2(m; — M\) + I. Using the WDLF's determined in the present calculations for pure DA 
WD's and X WD =50%, at most N ~ 2 WD's should have been expected in the HDF (4.4 
arcmin 2 ) at V«I«28 for th = 14 Gyr, whereas less than 1 is expected with the FORS1 
VLT telescope (V = 26; 6.8 arcmin 2 ). In the total available field of the french EROSII 
survey (250° 2 ), about iV ~ 26, 2 and 0.01 WD's are expected for halo ages th = 14, 15 
and 16 Gyr, respectively at Iu m = 20 in the appropriate (age-dependent) colors (see 
Fig. 2). These numbers are multiplied by a factor ~ 4 for I Hm = 21, a factor ~ 15 for 
Iii m = 22, and a factor ~ 1000 for J Km =25. Note that for th > 16 Gyr, the expected 
number of WD's remains essentially zero or a few. A useful tool for observers is the 
so-called discovery function D(M X ), i.e. the number of WD's observable over the whole 
sky. For a survey limited to nearby halo WD's, the density can be considered as constant, 
and D(M\) = ^-df nax (M\)n(Mx) . Figure 3 shows the expected discovery functions in the 
V-band for Xwd = 0.5 and different halo ages and magnitude limits. 

5. Conclusion 

We have shown in this paper that most of identified baryonic components are unlikely 
to provide a substantial contribution to the Galactic missing mass, except if they are 
distributed inhomogeneously, which implies a varying mass-to-light ratio in the Galactic 
dark halo. White dwarfs, although raising important issues about the Galaxy formation 
and evolution, remain the most plausible candidates to explain the observed microlensing 
events and might provide a substantial fraction of the sought baryonic dark matter. The 
luminosity functions, discovery functions and star counts calculated in the present paper 
will guide ongoing and future observational projects at faint magnitudes. As shown in Fig. 
2, the photometric observation of halo WD's will provide a powerful diagnostic to determine 
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the age of the Galactic halo. If WD's do indeed account for a large fraction of the Galactic 
missing mass, they solve the dilemna of the missing baryons at z — 0. In that case the 
baryonic missing mass is composed essentially of these stellar remnants in galactic halos 
and of the intergalactic left-over gas from the progenitors. 
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FIGURE CAPTIONS 

Fig. 1. — Expected halo WDLF for pure H-atmosphere WD's in the V-band, for three 
different halo ages and for X WD =50%. Solid line : IMF1; dashed line : IMF2 for t h = 14 
and 16 Gyr. The solid circles correspond to the faint end of the disk WDLF (Liebert et al., 
1998) while the squares correspond to the WD's in the sample with v tan > 250 km. s -1 

Fig. 2. — My vs (V-I) color- magnitude diagram for three different WD masses scaning the 
standard C/O WD mass range with ages in Gyr indicated by the filled circles for the mean 
value 0.8 M 

Fig. 3. — Discovery functions in the V-band for 2 magnitude limits, namely my = 22 (top 
lines) and m v = 20 (bottom lines), and for 2 halo ages, for X WD =50%. Solid and long-dash 
lines : IMF1; dot and short-dash lines : IMF2. Note that a kinematically selected sample 
might reduce the explored volume and thus the density of observable WD's as a function of 
the magnitude. 



